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We construct a model for inelastic neutrino- and electron-nucleon scattering cross sections using effective leading 
order parton distribution functions with a new scaling variable £ w . Non-perturbative effects are well described 
using the £ w scaling variable, in combination with multiplicative K factors at low Q 2 . Our model describes 
all inelastic charged lepton-nucleon scattering (including resonance) data (HERA/NMC/BCDMS/SLAC/JLab) 
ranging from very high Q 2 to very low Q 2 and down to the photo-production region. The model describes existing 
inelastic neutrino-nucleon scattering measurements, and has been developed to be used in analysis of neutrino 
oscillation experiments in the few GeV region. 



1. Introduction 

The field of neutrino oscillation physics has 
progressed from the discovery of neutrino oscil- 
lation [1 to an era of precision measurements of 
mass splitting and mixing angles. The cross sec- 
tions for neutrino interactions in the few GeV re- 
gion are not well known. This results in large 
systematic uncertainties in the extraction of mass 
splitting and mixing parameters in neutrino oscil- 
lations experiments such as MINOS 2J, NO^A[3] , 
K2K 4 , SuperKE], T2K0, and MiniBooNE[7]. 
A reliable model of neutrino inelastic cross sec- 
tions at low energies is essential for precise neu- 
trino oscillations experiments. 

The renewed interest in neutrino interactions 
at low energies has resulted in the construction 
of several near detectors (e.g. MINOS0, T2K[5]) 
to measure low energy cross sections and fluxes. 
In addition, there are new experiments (e.g. Sci- 
BooNE 0, and MINERvA[5]) which are specifi- 
cally designed to measure neutrino cross sections 
at low energies. 

In this communication, we report on a dual- 
ity based model of neutrino interactions. Ear- 
lier versions of the model (TOl 11] have been in- 
corporated into several Monte Carlo genera- 



tors of neutrino interactions including NEUT[12 , 
GENIE[13], NEUGENpJ] and NUANCE[H]. 

In the few GeV region, the three types of neu- 
trino interaction processes include quasi-elastic 
reactions, resonance production, and deep inelas- 
tic scattering. It is quite challenging to disen- 
tangle each of those contributions separately, and 
in particular the contribution of resonance pro- 
duction and inelastic scattering continuum. At 
low Q 2 there are large non-perturbative contri- 
butions to the inelastic cross section. These in- 
clude kinematic target mass corrections, dynamic 
higher twist effects, and higher order Quantum 
Chromodynamic (QCD) terms. 

In our previous studies [16117118] . non- 
perturbative effects were investigated within 
Leading Order (LO), Next-to-Leading Order 
(NLO) and Next-to-Next Leading Order (NNLO) 
QCD using charged lepton-nucleon scattering ex- 
perimental data [19 20121] , We found that in 
NLO QCD, most of the empirical higher- twist 
terms needed to obtain good agreement with the 
low energy data for Q 2 > 1 (GeV/c) 2 origi- 
nate primarily from target mass effects and the 
missing NNLO terms (i.e. not from interactions 
with spectator quarks). If such is the case, then 
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these terms should be the same in charged lep- 
tons (e, /i) and neutrino (i/^) scattering. There- 
fore, low energy data can be described by 
effective Parton Distribution Functions (PDFs) 
which are fit to high Q 2 charged lepton-nucleon 
scattering data, but modified to include target 
mass and higher-twist corrections that are ex- 
tracted from low energy e//i scattering data. For 
Q 2 < 1 (GcV/c) 2 additional corrections for non- 
perturbative effects from spectator quarks are re- 
quired. These corrections can be parametrized 
as multiplicative K factors. For the vector part 
of the charged current interaction the K factor 
terms should be the same in and ej[i scatter- 
ing. 

Therefore, a model that describes electron and 
muon scattering can also be used to model neu- 
trino scattering. At high Q 2 the vector and ax- 
ial structure functions are expected to be the 
same. The only terms which are not constrained 
by muon and electron scattering are terms which 
describe the difference between axial and vector 
structure functions and form factors at low values 
of Q 2 



2. Electro n-nucleon 
scattering 



and muon-nucleon 



In this section we define the kinematic variables 
for the case of charged lepton scattering from neu- 
trons and protons. The differential cross section 
for scattering of an unpolarized charged lepton 
with an incident energy E , final energy E' and 
scattering angle 9 can be written in terms of the 
structure functions T\ and T 2 as: 



(To- 
dQ.dE 



j(Eo,E', 9) — 



_ 4a 2 E'- 
Q 4 



cos 2 (6»/2) 



x [F 2 {x, Q 2 )/v + 2 tan 2 (0/2) ^(x, Q 2 )/M] 

where a is the fine structure constant, M is 
the nucleon mass, v = E Q — E' is energy of 
the virtual photon which mediates the interac- 
tion, Q 2 = 4E E' sin 2 (#/2) is the invariant four- 
momentum transfer squared, and x — Q 2 j2Mv is 
a measure of the longitudinal momentum carried 
by the struck partons. 

Alternatively, one could view this scattering 
process as virtual photon absorption. Unlike the 



real photon, the virtual photon can have two 
modes of polarization. In terms of the cross sec- 
tion for the absorption of transverse (or) and 
longitudinal (cr^) virtual photons, the differential 
cross section can be written as, 



d 2 a 
dQ.dE' 

where 



K 



7 =T[a T (x,Q 2 ) + ea L (x,Q 2 )] 
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The quantities T and e represent the flux and 
the degree of longitudinal polarization of the vir- 
tual photons respectively. The quantity R, is de- 
fined as the ratio ctl/ctt, and is related to the 
structure functions by, 



u( n2\ °l F 2 4M 2 x 2 T L 
R{x,Q z ) = — = ——{1+ )-l = — — 

a T 2xTx Q z 2xJi 



(5) 



where Fl is called the longitudinal structure 
function. The structure functions are expressed 
in terms of ol and ot as follows: 



MK 



Ti = 



An 2 a 

vK(cr L + <T T ) 

47r 2 a(l 



Q- 



T L {x,Q 2 ) = 1 + 



4M 2 i 2 - 

4M 2 x 2 



Q 2 



(6) 
(7) 

2xFi (8) 



or, 



2xF 1 =F 2 [l + ^^\-F L {x,Q 2 ). (9) 

Standard PDFs are extracted from global fits 
to various sets of deep inelastic (DIS) scattering 
data at high energies and high Q 2 , where non- 
perturbative QCD effects are negligible. PDF fits 
are performed within the framework of QCD in 
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either LO, NLO or NNLO. Here, using a new scal- 
ing variable (£ w ) we construct effective LO PDFs 
that account for the contributions from target 
mass corrections, non-perturbative QCD effects, 
and higher order QCD terms. 



The basic model: 
GRV98 PDFs . 



First iteration with 



Our proposed scaling variable, t; w is derived as 
follows. Using energy momentum conservation, 
the factional momentum, £ carried by a quark in 
a proton target of mass M is 



where 



2xQ 



'2 



Q 2 (l + v /l + 4M 2 x 2 /Q 2 )' 



(10) 



2Q' 2 = [Q 2 + Mf 2 ~ M 2 ] 

+ \J(Q 2 + M f 2 - Af, 2 ) 2 + 4Q 2 (M,; 2 + P 2 ) 

Here Mj is the initial quark mass with average 
initial transverse momentum Pr, and Mf is the 
mass of the final state quark. This expression 
for £ was previously derived [53] for the case of 
P T = 0. 

Assuming Mi = we construct following scal- 
ing variable 



2x(Q 2 + M f 2 + B) 



Q 2 [l + + AM 2 x 2 IQ 2 \ + 2 Ax ' 
or alternatively 

(Q 2 + M f 2 + B) 



Mv[l+ ^/l + Q 2 /v 2 ]+A' 



(11) 



(12) 



where in general Mf — 0, except for the case 
of charm-production in neutrino scattering for 
which Mf = 1.32 {GeV/c) 2 . The parameter A 
is used to account (on average) for the higher or- 
der QCD terms and dynamic higher twist in the 
form of an enhanced target mass term (the effects 
of the proton target mass is already taken into ac- 
count in the denominator of £ w ). The parameter 
B is used to account (on average) for the initial 
state quark transverse momentum, and also for 
the effective mass of the final state quark origi- 
nating from multi-gluon emission. A non-zero B 



also allows us to describe data in the photopro- 
duction limit (all the way down to Q 2 =0). 

In leading order QCD (e.g. GRV98 PDFs), 
J~2,lo for the scattering of electrons and muons 
on proton (or neutron) targets is given by the 
sum of quark and anti-quark distributions (each 
weighted the square of the quark charges): 



2,LO 



(x, t 



E^ 2 [xqi(x, Q 2 ) + xq { (x, Q 2 )] . (13) 



Our proposed effective LO PDFs model includes 
the following: 

1. The GRV98 LO Parton Distribution Func- 
tions (PDFs) ;23j are used to describe 
.Fg lo( x i Q 2 )- The minimum Q 2 value for 
these PDFs is 0.8 (GeV/c) 2 . 

2. The scaling variable x is replaced with the 
scaling variable £ w as defined in Eq. [TT| 
Here, 



3~ 2,Lo(* 2 '' Q ) ~ ^i e i 



(14) 



3. As done in earlier non-QCD based fits [2~4"l 
I26I27| to low energy charged lepton scatter- 
ing data, we multiply all PDFs by vector 
K factors such that they have the correct 
form in the low Q 2 photo-production limit. 
Here we use different forms for the sea and 
valence quarks, separately; 



K sea (Q ) 



Q 2 



j^vector (s~\2\ 



^valence 



(Q 



Q 2 + C s 

[1-G 2 D (Q 2 )] 
Q 2 + C„2 
Q 2 + C„i 



(15) 



where G D = 1/(1 + Q 2 /0.71) 2 is the proton 
elastic form factor. This form for the K 
factor for valence quarks is motivated by the 
closure arguments [28] and the Adler [29] 
US sum rule. At low Q 2 , [1 - G 2 D {Q 2 )} is 
approximately Q 2 /(Q 2 + 0.178). which is 
close to our earlier fit result [TU] . 
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These modifications are included in order 
to describe low Q 2 data in the photopro- 
duction limit (Q 2 =0), where J 7 ^ (x,Q 2 ) is 
related to the photoproduction cross section 
according to 



0.112 mb 



Q 2 



(16) 



A 

0.419 



B 

0.223 



0.554 



~C^2 X 2 /ndf 
0.31 1235/1200 



N 



valence 



1.011 [1-G 2 D (Q 2 



C 1 

w sea 

0.380 
Table 1 

First Iteration with GRV98 PDFs: Vector Pa- 
rameters. Only inelastic electron and muon scat- 
tering on hydrogen and deuterium (in the contin- 
uum region W > 2 GeV/c 2 ) are used in the fit 
Here the parameters are in units of (GeV/c) 2 . 



4. We freeze the evolution of the GRV98 PDFs 
at a value of Q 2 = 0.80 (GeV/c) 2 . Below 
this Q 2 , Ti is given by; 

T^{x,Q 2 < 0.8) = 

l^vector / s~)2\ -pvalence f £ r^l A Q\ 

-t^valenceW ) J ~2,LO \$v>iW — U.OJ 

+K vecto r(Q 2 ):F sea o ^ wjQ 2 = Qg) 

5. Finally, we fit for the parameters of the 
modified effective GRV98 LO PDFs (e.g. 
£ w ) to inelastic charged lepton scattering 
data on hydrogen and deuterium targets 
(SLAC PU /BCDMS [20] /NMC [2T]/H1[31]. 
In this first iteration, only data with an 
invariant final state mass W > 2 GeV/c 2 
are included, where W 2 = M 2 + 2M v-Q 2 . 

We obtain an excellent fit with the follow- 
ing initial parameters: ^4=0. 419, £?=0.223, and 
C„i=0.544, C„ 2 =0.431, and C sea =0.380, with 
X 2 /DOF = 1235/1200. Because of these addi- 
tional K factors, we find that the GRV98 PDFs 
need to be scaled up by a normalization factor 
iV=1.011. Here the parameters are in units of 
(GeV/c) 2 . These parameters are summerised in 
Table [U 

Thus, in our first iteration we modify the 
GRV98 T2 to describe low energy data down to 
photo-production limit as follows: 

F^(x,Q 2 ) = Q2 ® 38Q JITo(^,Q 2 ) 

+(1 ~ ^) Q2+o!544 ^" Ce( ^ lQ2)l (18) 
where £ w = . ^(Q 2 +Q-223) _ 

Q 2 ll + y/l + (2Mx) 2 /Q 2 ]+2*0.419x 



In fitting for the effective LO PDFs, the struc- 
ture functions data are corrected for the rela- 
tive normalizations between the SLAC, BCDMS, 
NMC and HI data (which are allowed to float 
within the quoted normalization errors). A sys- 
tematic error shift is applied to the BCDMS data 
to account for the uncertainty in their magnetic 
field, as described in the BCDMS publication|20). 
All deuterium data are corrected with a small cor- 
rection for nuclear binding effects [16|17|18] as de- 
scribed in section 12. We also include a separate 
additional charm production contribution using 
the photon-gluon fusion model in order to fit the 
very high energy HERA data. This contribution 
is not necessary for any of the low energy compar- 
isons, but is necessary to describe the very high 
energy low Q 2 HERA T2 and photoproduction 
data. The charm contribution must be added sep- 
arately because the GRV98 PDFs do not include 
a charm sea. Alternatively, one may use a charm 
sea parametrization from another PDF. 

The first iteration fit successfully describes all 
inelastic electron and muon scattering data in 
the continuum region including the very high 
and very low Q 2 regions. We find that al- 
though photo-production data were not included 
in our first iteration fit, the predictions of our 
model for the photo-production cross sections 
on protons and deuterons (Q 2 = limit) are 
also in good agreement with photoproduction 
measurements [32] . 

Furthermore, although no resonance data were 
included in the first iteration fit, the fit also pro- 
vides a reasonable description of the average value 
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of J-~2 for SLAC and Jefferson data in the reso- 
nance region [33] (down to Q 2 = 0.07 (GeV/c) 2 ). 

4. Second iteration with GRV98: Includ- 
ing photo-production data, resonances, 
and additional parameters 

We now describe the second iteration of the 
fit [IT]. Theoretically, the K t factors in Eq. 15 
are not required to be the same for the u and d 
valence quarks or for the u, d, s, sea quarks and 
antiquarks. In order to allow flexibility in our 
effective LO model, we treat the Ki factors for u 
and d valence and for sea quarks and antiquarks 
separately. 

In this second iteration, in order to get ad- 
ditional constraints on the different Ki factors 
for up and down quarks separately, we include 
photo-production data in the continuum region 
(W > 2 (GeV/c) 2 ) for both hydrogen and deu- 
terium. We also include electron scattering data 
in the resonance region (on hydrogen and deu- 
terium) in the fit. In order to extract neutron 
cross section from photproduction cross sections 
on deuterium, we apply a small shadowing correc- 
tion as shown in figure [T] The small nuclear bind- 
ing corrections for the inelastic lepton scattering 
data on deuterium is described in section \Wl 



LW 
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sea — strange ) 



p, _|_ (jL — Ehad 



ryi _|_ /^vector 

* sea— strange 

Q 2 

r^2 _i_ (^vector 
^fr ' ^ sea— up 

Q 2 

Q2 _|_ (^vector 
sea— down 

KvlllZce-upiQ 2 ) = K LW [1 - G 2 D (Q 2 )} 

vecti 
v2u 



■^sea— up\ V ) 



T^vector //~)2\ 
K 



"■ sea—dowvi 



Q2 j_ /^vector 

q 2 + cm 



K 



vatence- 



-dowi 



(Q ) = K LW ([1-G 2 D (Q 2 )} 

(r\2 j_ /-^vector \ 
Q ) (1Q) 

Q 2 + c vtd or ) 

The best fit is given by A = 0.621 ±0.009, B = 
0.380 ± 0.004, C^ tor = 0.341 ± 0.007, C^ tor = 



1.02 



1.00 



0.98 



0.96 




0.2 0.5 1.0 2.0 5.0 10.0 20.0 
Pbeam (GeV) 



Figure 1. The ratio of photoproduction cross sec- 
tions on deuterium to the sum of the photopro- 
duction cross sections on unbound protons and 
neutrons. This shadowing correction is used to 
extract the photorproduction cross section on free 
neutrons and protons. 



0.417 ± 0.024, C™f or = 0.323 ± 0.051, C^ tor = 
0.264±0.015, and C L ~ Ehad = 0.217±0.015. The 
sea factors are ^^=0.621, C^ p =0.363, 
and C v s la° r strange was set to be the same as 
^sea-down ■ Here, the parameters are in units 
of (GeV/c) 2 . The fit x 2 /DOF =2357/1717, and 
N = 1.026±0.003. The resonance data add to the 
\ 2 /ndf because the fit only provides a smooth av- 
erage over the resonances. No neutrino data are 
included in the fit. These parameters are summa- 
rized in Tabled] 

The normalization of the various experiments 
were allowed to float within their errors with the 
normalization of the SLAC proton data set to 1.0. 
The fit yields normalization factors of 0.986 ± 
0.002, 0.979 ±0.003, 0.998 ±0.003, 1.008 ±0.003, 
1.001 ± 0.004, and 0.987 ± 0.005 for the SLAC 
deuterium data, BCDMS proton data, BCDMS 
deuterium data, NMC proton data, NMC deu- 
terium data, and HI proton data, respectively. 
With these normalization, the GRV PDFs with 
our modifications should be multiplied by N = 
1.026 ±0.003. 
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B 


C v 2d 


C v 2u 


0.621 


0.380 


0.323 


0.264 


^ sea 


^ sea 


G'vld 




0.561 


0.369 


0.341 


0.417 


/^strange 
^ sea 


(jL — Ehad 




&-A1232 
1/ — u 


0.561 


0.218 


1.65 


0.55 


X 2 /ndf 


*F valence 


N 




2357/1717 




1.026 





Table 2 

Second iteration with GRV98 PDFs: Vector Pa- 
rameters. Here, we now include electron scat- 
tering data in the resonance region and pho- 
toproduction data on hydrogen and deuterium. 
The resonance data adds to the x 2 /ndf because 
the fit provides a smooth average over the reso- 
nances. No neutrino data are included in the fit. 
When applicable, all parameters are in units of 
(GeV/cf. 



Note that we apply also a small d/u correction 
to the GRV98 PDFs as described in section HU 
This correction increases the d quark distribution 
at large x and is extracted from NMC data for 

As seen in the plots, our fit describes all of 
the data, including photoproduction data in the 
continuum region. In the resonance region, the 
fit works well, except near Q 2 = 0. In or- 
der to better describe the resonance region near 
Q 2 = we include the low W K LW factor, where 

K LW = ( v 2 + C L-Ehady( v 2y 

For neutrino scattering we must treat the I — 
3/2 A1232 resonance as a special case for all 
values of Q 2 . This is because the neutrino- 
proton charged current scattering amplitude for 
the production of the A(1232) / = 3/2 resonance 
is y3 times the amplitude in neutrino-neutron 
scattering 34 35J . Similarly, the antineutrino- 
neutron charged current scattering amplitude for 
production of the A (1232) 7 = 3/2 resonance 
is v3 times the amplitude for the production 
in neutrino-proton scattering. This relationship 
holds for both the vector and axial amplitudes. 

Therefore, for the case of charged current neu- 
trino and antineutrino scattering the K LW fac- 
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Figure 2. The effective LO PDF model com- 
pared to charged lepton Ti experimental data 
(SLAC, BCDMS, NMC) at high x (these data 
are included in our fit) :[top] T2 proton, [bot] 
T<x deuteron. The solid lines are our fit, and the 
dashed lines are GRV98 . 



tors for the valence u and d quark distributions 
are modified as follows: 



K, 



LW 



C 



L — Ehad 



(W > 1.4 GeV/c 2 ) 



k lw 

vu 
K LW 

j^LW 



C 



L — Ehad 



= K, 
= K, 



v 

A1232 



A1232 



(W > 1.4 GeV/c 2 ) 

= 1.65 (W < 1.4 GeV/c 2 ) 
= 0.55 (W < 1.4 GeV/c 2 pO) 



Comparisons of our fit to various sets of in- 
elastic electron and muon Ti data on proton 
and deuteron targets are shown in Figures |2j (for 
SLAC, BCDMS and NMC). Comparisons to low 
x HI (electron) data on protons are shown in Fig- 
ure [3] Our effective LO model describes the in- 
elastic charged lcpton F 2 data both in the low x as 
well as in the high x regions. The model also pro- 
vides a good description of both low energy and 
high energy photo-production cross sections 32J 
on proton and deuteron targets as shown in Fig- 
ure [SI 

5. Comparison to resonance production 
data 

Comparisons of the model fit to hydrogen and 
deuterium electron scattering data in the reso- 
nance region |33j are shown in Figure |4j As ex- 
pected from quark-hadron duality [36] , our model 
provides a reasonable description of both the in- 
elastic region as well as the average value of 
the T% data in the resonance region (down to 
Q 2 = 0.0), including the region of the first res- 
onance (W = 1.23 GeV/c 2 )). We find also good 
agreement with the most recent !Fl and T<x data 
in the resonance region from the E94-110, and 
JUPITER experiments [33|37j at Jlab, as shown 
in Fig. [5] Our predictions for Fl are obtained us- 
ing our T<x model and the -ftiggs [38] parametriza- 
tion (as discussed in section 9). We find good 
agreement with quark hadron duality down to 
very low Q 2 . Other studies [30 with unmodi- 
fied GRV PDFs find large deviations from quark- 
hadron duality in the resonance region for elec- 
tron and muon scattering. This is because those 
studies do not include any low Q 2 K factors and 
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Figure 3. The effective LO PDF model compared 
to charged lepton Ti experimental data at low x 
from HI (these data are included in our fit). 



use the scaling variable £ (while we use the mod- 
ified scaling variable £ w ). We find that quark 
hadron duality works at low Q 2 if we use the mod- 
ified scaling variable £ w , and low Q 2 Ki factors. 

In the Q 2 = photoproduction limit, the 
model provides a good descriptions of the data 
for both the inelastic region as well as in the res- 
onance region as shown in Figure [6] Since the 
photoproduction cross sections on protons and 
neutrons is are not far aprt, our fits yield the 
interesting result that at Q 2 — 0, the effective 
valence d quark is multiplied by a factor of 0.85 
(=0.323/0.341,) and the effective valence u quark 
distribution at Q 2 — 0) is multiplied by a fac- 
tor of 0.634 (=0.264/0.417). This implies that at 
Q 2 = 0, the effective d/u ratio for valence quark 
is enhanced by a factor of 1.5. 



8 



. 1 1 "|''"I'"'I"''I''TTT'I 
- Cnu=0.21B 




— - Csnu=0 




Q 2 =0.07 


i 






I....I.-.I--I..-I.-I-.U 





0.02 0.03 0.05 0.07 0.10 
X 





0.0010 
0.0005 





Q 2 =9.0- 








V 







0.4 
0.3 
0.2 
0.1 


0.15 
0.1 
0.05 

J o 

0.1 
0.05 


0.075 
0.05 
0.025 


-0.025 



Q =1.5 GeV 2 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



E=2.4 GeV 

theta=20 

0.2<Q 2 <0.6 



0.2 0.4 
X 





................. I ... . 

E=4.0 GeV 




theta=30 




1.6<Q 2 <2.7 


— I . . 


. . i — i — i — i 



0.3 0.4 0.5 0.6 0.7 0.3 0.9 
X 







E=3.2 GeV 






theta=27 






k 0.8<Q 2 <1.6 










. . 1 






Figure 5. Comparisons of the predictions of our 
model to proton data for Tl (note that data for 
Th are not included in our fit). 



6. Application to neutrino scattering 

For high energy neutrino scattering on quarks 
and antiquarks, the vector and axial contribu- 
tions are the same. At very high Q 2 , where the 
quark parton model is valid, both the vector and 
axial K factors are expected to be 1.0. There- 
fore, high Q 2 neutrinos and antineutrino struc- 
ture functions are given by : 



Figure 4. Comparisons of charged lepton exper- 
imental data in the resonance region to the pre- 
dictions of our effective LO model: [top] six pro- 
ton plots, [bot] four deuteron pots. The red line 
includes the K LW factor and the dashed black 
line does not include the K LW factor. Note that 
the K LW factor is not used in the region of the 
A(1232) resonance. 



^(x,Q 2 ) = 2X l [t w q l (Z w ,Q 2 )+Z w q l (Z Wl Q 2 )]. 
and 

xT^(x,Q 2 ) = 2J: i ^ w q^ w ,Q 2 )~^ w q l (iw,Q 2 )]- 
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Figure 6. The effective LO PDF model com- 
pared to photoproduction cross sections (Q 2 = 
limit) at low and high energies (these data are 
included in our fit), [top] proton, [bot] deuteron. 
At very high photon energy, we include charm 
contribution from gluon fusion process which is 
needed to describe the very high energy HERA 
data. If we want to also describe the photopro- 
duction data in the resonance region, we need to 
multiply the u and d valence PDFs by K LW = 
(v 2 + C L - Ehad /{v 2 + 2Q 2 ). The red line includes 
include the K LW factor, and the dashed black 
line does not includes the K LW factor. Note that 
the K LW factor is not used in the region of the 
A(1232) resonance. 



where 








q up = 


d + s; 




— u + c 


q vn = 


u + s; 


q vp 


= d + c 


q Vp = 


u + c; 


q w 


= d + s 


q Vn = 


d + c; 


q vp 


= u + s 



(21) 

Here, T 2 = t>W 2 , T\ = MWi and J" 3 = ^W 3 . 
Note that for the strangeness conserving (sc) part 
of the u and d quark distributions, the PDFs 
are multiplied by a factor of cos 2 9 c . For the 
strangeness non-conserving part the PDFs are are 
multiplied by a factor of sin 2 9 c . 

The neutrino structure functions must satisfy 
the following inequalities: 




Wi 



< 



> 



Vo 2 



Vq 2 



72 \m 

2M 

72 \m 

2M 



Wi < (1 + v 2 /Q 2 )W 2 



(22) 



There are several major difference between the 
case of charged lepton inelastic scattering and 
the case of neutrino scattering. In the neutrino 
case we have one additional structure functions 
J-^(x,Q 2 ). In addition, at low Q 2 there should 
be a difference between the vector and axial Ki 
factors due a difference in the non-perturbative 
axial vector contributions. Unlike the vector T 2 
which must go to zero in the Q 2 = limit, we ex- 
pect |24|40j that the axial part of J- 2 is non-zero 
in the Q 2 = limit. 

We already account for kinematic and dynamic 
higher twist and higher order QCD effects in 
T 2 by fitting the parameters of the scaling vari- 
able £ w and the K factors to low Q 2 data for 
J-2n(x,Q 2 ). These should also be valid for the 
vector part of T 2 in neutrino scattering. How- 
ever, the higher order QCD effects in T% and xJ-3 
are different. We account for the different scaling 
violations in J- 2 and by adding a correction 
factor H(x, Q 2 ). 

The axial K factors and the H correction to 
xT-i are included in the following expressions: 



"(as,' 



f) = ^Kr ctor (Q 2 )^q^ w ,Q 2 ) 

' tor (Q 2 )^U^,Q 2 )m) 
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Figure 7. The x and Q 2 dependence of the factor 
H(x, Q 2 ) that accounts for the difference in the 
QCD higher order corrections in T<x and xT 3 



and 

^ axial (x,Q 2 ) = ZiKr al (Q 2 )£ w qi(t w ,Q 2 ) 

+X j K? xial (Q 2 )Z w q j (Z w ,Q 2 ) (24) 

and 

xT^(x,Q 2 ) = 2H(x, Q 2 ) J ^iKf^ wqi ^ w ,Q 2 ) 



-ZjKfz^qjiUQ 2 )} (25) 



where, i denotes (valence — up), (valence 
down), (sea — up), (sea — down), and (sea 
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Figure 8. A fit to the x dependence of the factor 
H(x, Q 2 ) that accounts for the difference in the 
QCD higher order corrections in Ti and (at 
Q 2 =8 (GeV/c) 2 ). 



strange), and j denotes (sea — up), (sea — down), 
and (sea — strange) . The factor Kf F3 is given by 



K 



xF3 



Q2K? xml {Q 2 )Kf ector (Q 2 ) 



C 



xF3 



(26) 



The parameter C x fz = 0.33 is set such that the 
Adler sum rule for xT^ is satisfied at Q 2 = 0. 

We obtain an approximate expression for 
H(x, Q 2 ) as the ratio of two ratios as follows: 
H(x,Q 2 )= D xF3 (x,Q 2 )/D F2 (x,Q 2 ) where 



D xF3 (x,Q 2 ) 
D F2 (x,Q 2 ) 



xTS l °(x,Q 2 ) 
xT l 3 °(x,Q 2 ) 

T^ lo (x,Q 2 ) 
T\°(x,Q 2 ) 



(27) 



The double ratio H(x,Q 2 ) is calculated by the 
TR-VFS scheme[39] with MRST991 NLO PDFs. 
This ratio turns out to be almost independent 
of Q 2 . The results of this calculation at Q 2 = 
8(GeV/c) 2 are fit by the following functional 
form: 



H(x,Q 2 ) = 0.914 + 0.296x 

- 0.374x 2 + 0.165a; 3 



(28) 
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We use the above function as an approximation 
for H(x,Q 2 ) for all values of Q 2 . 

In our previous |10lll| analysis we assumed 
H(x,Q 2 )=l, C xF3 = 0, and Kf xial {Q 2 ) = 
K^ ector (Q 2 ) . However, this assumption is only 
valid for Q 2 > 1 (GeV/c) 2 . Here, we improve on 
our previous analysis by introducing Kf x (Q 2 ) 
factors which are different from Kf ector (Q 2 ) ) and 
include the H(x,Q 2 ) and C x f 3 corrections for 
xT 3 . 



jy-axial 

sea— strange 



K 



sea~-dow7i 



(Q 2 ) 



(Q 2 ) - 



j^axial ( /~\2\ 



sea— up 



(Q) 



Q 2 
Q 2 + 

Q 2 
Q 2 + 



r> r^axial 

sea — strange 

I /~faxial 

' ^ sea— strange 

D r^axial 
« sea— down 
axial 



Q 2 



r> r^axial 
J-u^ sea— up 



r^axial 
^ sea— up 



In addition, as mentioned earlier 



(30) 



7. Axial structure function T<x, and the 
structure function xT 3 

The axial structure function T% is not con- 
strained to go to zero at Q 2 = 0. At higher Q 2 
(e.g. Q 2 > 1) the vector and axial structure func- 
tions should be equal. 

The simplest assumption to make is that for 
valence quarks the parameters of parton distri- 
butions and associated K factors for the axial 
case are the same as the parameters for the vec- 
tor case, except for the parts that constrain Ti to 
be zero at Q 2 . Here, we assume that the axial K 
factor for valence quarks is equal to 1 (versus the 
vector valence K factor which is (1 — Gd(Q 2 ))- 

Therefore, the suggested axial K factors are: 



K. 



LW 



C 



L-Ehad 



LW 



K. 



K LW 
LW 



C L 



(W > 1.4 GeV/c 2 



•Ehad 



K. 



and 



v 

K v-d 
j^A1232 



, — (W > 1.4 GeV/c 2 ) 

-- 1.65 {W < 1.4 GeV/c 2 ) 
-- 0.55 (W < 1.4 GeV/c 2 j29) 



K 



axial 

valence— dow 

ly-axial i (~)2\ 

valence— up V^v ) 



n(Q 2 ) = K, 



-LW 1 


f Q 2 - 


i_ (^vector 
h ^v2d 


vd \ 


\Q 2 - 


L /^vector 
h ^vld 




( Q 2 - 


l_ /^vector 
r ^v2u 


-LW | 




vu \ 


kQ 2 - 


i_ s~ivector 
l" ^vlu 



For sea quarks, we need to include additional pa- 
rameters 



K: 



xF3 



Q 2 Kf xlal {Q 2 )K^ ector {Q 2 ) 
Q 2 + C xF3 



(31) 



where the parameter C X F3 = 0.33. With this 
value, we get that W 3 (Q 2 = 0) for the A (1232) 
is equal to 11. This is close to the theoretical 
calculation of of Lalakulich and Paschos[34 . In 
addition, with this value the Adler rule for xT 3 
is approximately satisfied at Q 2 = 0. 
We make the following assumptions 



P = P u =Pd = P S = 

s~iaxial 

^ sea— strange 

s~iaxial 

sea— down 



C: 



axial 
sea— up 



0.6 ±0.3 
0.3 

= 0.3 
= 0.3 



(32) 



With these choices the axial K factors are be- 
tween K = 1 (P=l) and K axial = K vector . Here 

P = 0.6 ± 0.3 yields J^ xial \^ w = 0.00001, Q 2 = 
0)( p+n )/2 = 0.33 ± 0.16 for the average of the 
neutron and proton structure functions, which 
is the value calculated from PCAC in the model 
of Kulagin and Peti|40j. For an iron target (as- 
suming a ratio Fe/D = 0.8) this corresponds to 
j-axtal)^ = 0.00001, Q 2 = 0) Fe = 0.25 ±0.11, 
which is similar with the value of 0.210 ± 0.02 
measured by CCFR [25] using a different func- 
tional form for the extrapolation of their data to 
to Q 2 = 0. The value of C^ al = 0.3 was se- 
lected because it yields a PCAC contribution for 
j-axtai)^ = o ooooi, Q 2 = 1) = 0.08, which is 
the value calculated in the model of Kulagin and 
PetigO]. 

A better determination of the axial K factors 
for the sea quarks requires including neutrino 
data on nuclear targets in our fit. The nuclear 
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corrections at small x for vector and axial struc- 
ture functions are different in general. Therefore, 
it is not possible to separately determining the 
axial K factor without assuming a specific model 
for the neutrino nuclear corrections 40 . There- 
fore, this study will be presented in a future com- 
munication. 



8. Comparison to neutrino and antineu- 
trino data on the production of the A 
1232 Resonance 

Our model is aimed at describing neutrino 
cross sections in the inelastic region (e.g. 
W > 1.9 GeV/c 2 ) and match on to other 
calculations (34135] that describe the resonance re- 
gion and the quasielastic peak. In order to match 
smoothly to a resonance region model, the model 
should also describe the average of the cross sec- 
tion in the resonance region. Since the model 
describes electron and muon scattering structure 
functions in the resonance region (on average, as 
expected from local duality), it should also work 
for the vector part of the structure functions in 
neutrino scattering. 

We now test how well the suggested axial pa- 
rameters for the valence quarks describe the pro- 
duction of the A(1232) resonances on proton and 
neutron targets. 

At a fixed value of the final state invariant mass 
W, the differential cross section for neutrino scat- 
tering at incident energy E is given [3"0] by: 



da 



G 2 



dQ 2 dW 4tt 



cos tic 



' M 2 



w 

ME 2 
1 



W^ + m 2 ) 



2(k-p)(k> ■p)-^M 2 (Q 2 +ml) 



(33) 



+ 



M 2 



Q 2 k-p - \q-p{Q 



2 i 2 \ 



W 4 . (Q 2 + ml) 



M 2 » 



5 2 



M 2 



mJk-p) 



where k and k' are the momenta of the initial and 
final state lepton, and p is the momentum of the 



target nucleon. This expression can be written as 

da G 2 , „ W I 1 

cos Be 



dQ 2 dW 4tt 



-Wo 



Wi(Q 2 + m^) 



2(1- 



m 2 

Jil 2 




vector part of W4 and W5 are related to the 
vector part of W2 and W\ by the following 
expressions 34 : 



yyvector -yyvector 



Mv 
M 2 v 2 



Wi 



M 2 

Q 2 



(35) 



The contribution of the W4 and W5 should 
be accounted for in the case of r neutrino 
interactions [41], For muon neutrinos, we find that 
even for a neutrino energy as low as 1 GeV, the 
contribution of the W4 term is negligible. Us- 
ing the model of O. Lalakulich, E. A. Paschos 
for the production of the A (1232) resonance[34 
we find that for a neutrino energy of 1 GeV, the 
contribution of the W5 term is about 10% for Q 2 
close to m^. We also find that the assumption 

yyaxial _ yyvector ^ and yyaxtal = yyvector jg good 

to within 30% for Q 2 = 0. 

At high energy (E > 20 GeV), for low values 
of W (for which v is small) and low Q 2 (Q 2 < 
3 (GeV/c) 2 ), all the terms are small except for W2 
term. For this case the above expression reduces 
to : 

• " -(W,Q 2 ) 



dQ 2 dW 



G 2 2n WT% N ^ A 

— cos^ B c — 

2tt M 



cos M v 



(W, 



(36) 



where S cos = ^- cos 2 B c = 80 x 10 40 cm 2 , and 



N A is the average value of F2 in the region of 
the A (1232) resonance (1.1 < W < 1.4 GeV/c 2 ). 
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Figure 9. 

A (Q 2 ) 



'nun 

•'a 



12 3 

Q 2 (GeV/c) 2 

A comparison of J 7 ; 
and the sum T', 



v(p+n)- 



vp— >A 
2 

A,^2 



(Q 2 ), 



\Q 2 ) 



tracted from A production data to our predic- 
tions. Solid line, our updated model with sug- 
gested axial parameters, and K A } 232 = 0.55 and 
#A1232 = ]_ 6g The dashed blue Une is predic- 
tions of our model with the assumptions of Ax- 
ial= Vector and K A }f 2 = 1 and K A }f 2 = 1. 



For a specified range of the final state invariant 
mass Sw we get. 



-pvp- 



and 



•> 2 



l (Q 2 



\Q 2 ) 



vM da v P^ A 



ir 



dQ 2 



vM da 



s rn ,ws 



ir 



dQ 2 



(37) 



(38) 



For the region of the A (1232) resonance W — 
1.23 GeV/c 2 , and 6 W = 0.3 GeV/c 2 , which cor- 
responds to the invariant mass range from pion 
threshold M w +M = 1.1 GeV/c 2 to 1.4 GeV7c 2 . 

The antineutrino structure functions are re- 
lated to the neutrino structure functions by the 
following rations. 



-pvp 

•'A 

-pun 



-pvn 
— Tr v P 



(39) 



For da vp ^ A /dQ 2 in the 1.1 < W < 1.4 GeV/c 2 
region we take the average of da" p ^ p ' K+ /dQ 2 
data of Allen et al [32] and Allasia et al |4"3"] . 
Since T\ v = F» n we also include the All- 
sia et al data for da yn ^ n ' K ~ / dQ 2 in this aver- 
age. For da vn ^ A /dQ 2 data we sum the Allasia 
et al measurements for da vn ^ v ^ /dQ 2 and the 
dcr lm_m7r /dQ 2 measurements in the 1.1 < W < 
1.4 GeV/c 2 region. 

Figure [9] shows a comparison data for 
1.23, Q 2 ), F^ A {W = 1.23, Q 2 ) 



and for the sum F^ J ' +r ' 



1.23, Q 2 ), ex- 
tracted from the Allen et. al [32] and Allasia et 
al [43] measurements, to to our predictions (with 
K A } 23 2 = 0.55 and K A } 232 = 1.65. There is 
good agreement between the A production data 
and our model. 

Also shown is the predictions of our model 
under the assumption that Axial=Vector and 
K A }_ 232 = 1 and K A }f 2 = 1. Using the sug- 
gested axial parameters yields a better descrip- 
tion of the data near Q 2 = 0, and also a better 
description of the difference between the scatter- 
ing on protons and neutrons. 

For most applications the target consists of ap- 
proximately equal numbers of neutrons and pro- 
tons. For such isoscalar targets, our previous 
version which assumed that Axial=Vcctor and 
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extracted from the measured differential cross 
section for hyperon production from Fermilab 
experiment E180 [44] to our predictions with 
K LW (W < 1.4) = 1 . There is good agree- 
ment between the Ao, So production data and our 
model. 



^A_1232 = ! and ^A_1232 = j 



worked reason- 
ably well except near Q 2 = 0. Our model with 
updated axial parameters works much better. 



9. Comparison to quasielastic production 
of the Ao and So hyperon resonances. 

The expression for quasielastic production of 
Ao(1116) and So(1193) hyperon uds resonances 
in antincutrino-proton collisions is given by: 



dQ 2 dW 



g 2 ■ 2. wr 2 

— sin On 

2tt M 



up— yAo.Ho 



= s. 



W F. 

'■M~ 



i/p->A ,E 



(Q 2 



<Q 2 ) 

(40) 



where S sin = ~ sin 2 9 C 



4.29 x 10 



-40 



cm 



Here, Ti is the average value of for the range 
1.1 < W < 1.4 GeV/c 2 , where 1.4 GeV/c 2 is 
the proton-kaon threshold, and 8 W = 0.3 GeV/c 2 . 
Therefore, in the expressions are identical to the 
expressions for the production of the A1232, ex- 
cept that here we have sin 2 Qq instead of cos 2 6c- 
Since we do not have a / = 3/2 resonance, we use 
K LW (1.1 < W < 1.4) = 1. 

Figure [lO] shows a comparison of data for 
j~vp-tA ,a ex t rac ted from the measured differ- 
ential cross section for hyperon production from 
Fermilab experiment E180 44J to our predictions 
with K LW (W < 1.4) = 1 . There is good agree- 
ment between the hyperon (Ao,So) production 
data and our model. 



10. 



2xJ=l ector 
function 



and the longitudinal structure 



In the extraction of the original GRV98 LO 
PDFs, no separate longitudinal contribution was 
included. The quark distributions were directly 
fit to J-2 data. A full modeling of electron and 
muon cross section requires also a description of 
2x T\ . We use a non-zero longitudinal R in recon- 
structing 2xFi by using a fit of R to measured 
data. In general, IxT^ is given by 



2xT e ^{x,Q i ) = Fl'»(x,Q z )x 
1 + AM 2 x 2 /Q 2 



(42) 



l + R{x,Q 2 



or equivalent ly 



1 + v 2 jQ 2 



(43) 



l + R(x,Q 2 ) 



The i?i998 function [55] provides a good descrip- 
tion of the world's data for R in the Q 2 > 0.30 
(Gev/c) 2 and x > 0.05 region (where most of the 
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R data are available). 

R e / m u{x, Q 2 > 0.3) = #i99 8 (a:, Q 2 > 0.3) 



to the vector component of 2x JF\. Investigation 



(44) 



However, for electron and muon scattering (and 
for the vector part of neutrino scattering) the 
Ri 99 $ function breaks down below Q 2 — 0.3 
(Gev/c) 2 . Therefore, we freeze the function at 
Q 2 = 0.3 (Gev/c 2 ) and introduce additional pa- 
rameters for R in the Q 2 < 0.3 (Gev/c 2 ) region. 
For electron and muon scattering and for the vec- 
tor part of T\ we introduce a K factor for R in 
the Q 2 < 0.3 (Gev/c 2 ) region. The new function 
provides a smooth transition for for R e i il from 
Q 2 = 0.3 (Gev/c 2 ) down to Q 2 = by forcing 
Rvector to approach zero at Q 2 = as expected in 
the photoproduction limit. This procedure keeps 
al/Q 2 behavior at large Q 2 and matches to Ri 99 $ 
at Q 2 = 0.3 (Gev/c) 2 . 



R e/(I (i,O 2 < 0.3) = 3.633 x 



Q 2 



1 



x R 199S (x,Q 2 =0.3) 

Using the above fits to R as measured in elec- 
tron/muon scattering we use the following expres- 
sions for the vector part of 2a;J r i neutrino scatter- 
ing. 



2xTr ctor (x,Q 2 ) = jq 



"(a;, i 



f) x 



R 

vector 

(x,Q 2 > 0.3) 

Rvector(x, Q 0.3) 

or equivalently 



1 + 4M 2 x 2 /Q 2 
l + R(x,Q 2 ) 
R 199a (x,Q 2 > 0.3) 
Re/»(x,Q 2 <0.3) 



Wr c * or (x,Q 2 ) = W 2 " ector (x,Q 2 ) x (45) 

1 + v 2 /Q 2 

1 ~t~ ^vector (*^ ? Q ) 

The above expressions have the correct limit at 
Q 2 = 0. 

11. The structure function 2xJ r i xlal 

As is the case for F% = vW2, the axial com- 
ponent of 2xJ-\ (= 2xMW\) at large Q 2 is equal 



of the Adler sum rules (described in section 151 
shows that at Q 2 = 0, the Adler sum rules for the 
axial structure functions yy® xlal and y\l^ xtal are 
the same. 

Therefore, we assume that W? xial (Q 2 = 0) = 
W% xtal (Q 2 = 0). This also implies that the longi- 
tudinal part of y\?£ xzal dominates at Q 2 — as ex- 
pected from PCACgD]. Here, since T% xlal (Q 2 = 
0) is non zero, it implies that R axial (Q 2 = 0) is 
infinite. Therefore, it is more convenient to cal- 
culate Wf xlal instead of Ff xial at low Q 2 . 

In summary, at large Q 2 we expect: 

W? xial (x,Q 2 >1) = Wp ctor (x,Q 2 ) 

And at Q 2 = we expect: 

W?" al (x,Q 2 = 0) = W% xial (x,Q 2 ) 

We use a function for Y\;^ xlal that matches 
these two boundary conditions as follows. 



Wj (x,Q ) = W 1 (x,Q ) 



+ W% xial (x,Q 2 )(l- 



Q 2 
Q 2 + 0.2 

Q 2 



Q 2 + 0.2 



(46) 



12. Charm production in neutrino scat- 
tering 

Neutrino scattering is not as simple as the case 
of charged lepton scattering because of charm 
production. For the non-charm production com- 
ponents we use J-2 Cp (x,Q 2 ), 2xJ-" cp (x,Q 2 (sum 
of vector and axial parts) and xJ-^ cp (x, Q ) as 
described above. 

For the charm production components of 
T?\x,Q 2 ), xF^(x,Q 2 ) and 2xJ 7 ?\x,Q 2 ) the 
variable £ w now includes a non-zero M c = 
1.32GeV/c 2 . 

The target mass calculations as discussed by 
Barbieri et. al 22 imply that J 7 ^ ~ cp is described 
by J*2 ~ cp (^ w , Q 2 ), and the other two structure 
functions are multiplied by the factor K c h arm = 
q2+ M '2 ■ Therefore, to include charm production 
we use the following expression for charm produc- 
tion processes. 



1G 



0.001 



^■charrn 
2\ 



Q 2 



Ml 



Q 2 

? vector - cp (x,Q 2 ) = Z t Kr ctor (Q 2 

x [twqi{tw,Q 2 ) + £wqi(Zw,Q 2 )] 



■v, axial— cp 



(x, t 



Ix^F i {x.Q ) — K charm 



l+4M 2 x 2 /Q 2 
i + R(UQ 2 ) J - 



T?(x,< 



(47) 



(48) 



and 



x^(x,Q 2 ) = 2H(x,Q 2 )K charm 
-^Kf^ w q^ w ,Q 2 )\ 



(49) 



Here we use the -R1998 parametrization |19j for 
the vector part of R nc P and R cp 

13. Nuclear corrections 

In the comparison with neutrino charged- 
current differential cross section on iron, a nu- 
clear correction for iron targets should be applied 
applied. Previously, we used the following pa- 
rameterized function, fix) (a fit to experimen- 
tal electron and muon scattering data for the ra- 
tio of iron to deuterium cross sections, shown in 
Fig[TT|), to convert deuterium structure functions 
to (isoscalar) iron structure functions |45j : 

f(x) = 1.096 - 0.364 x- 0.278 e ~ 21 - 94 x 

+ 2.772 x 14 - 417 (50) 

However, we find that the ratios of iron to deu- 
terium structure functions from SLAC and Jef- 
ferson Lab are better described in terms of the 
variable £ w . In addition, if £ w is used, the func- 
tion that describes the iron to deuterium ratios 
in the deep inelastic region is also valid in the 
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Figure 11. The ratio of Ti data for 
heavy nuclear targets and deuterium as mea- 
sured in charged lepton scattering experi- 
ments(SLAC,NMC, E665). The band show the 
uncertainty of the parametrized curve from the 
statistical and systematic errors in the experi- 
mental data 1431. 



resonance region. Therefore, we use the following 
updated function F upd ated{tw)- 

Fupdatedi^w) = 1-096 - 0.364 £ w 
- 0.278e~ 21 - 94 «™ 



14.417 



(51) 



Figure [12] shows a comparison of Jefferson lab 
measurements of ratio of the cross sections on 
iron to deuterium in the resonance region[4"6"] to 
data from SLAC E87 and SLAC E139 in the deep 
inelastic region. The data are plotted versus x 
and are compared to the fit function fix"). Fig- 
ure [13] shows the same data plotted versus £ w 
and are compared to our updated fit function 
Fvpdated(€w)- When plotted versus £ w the ra- 
tios in the deep inelastic region and the ratios in 
the resonance region all lie on the same universal 
curve. 

For the ratio of deuterium cross sections to 
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a Jlab Resonance, ♦ SLAC E87, — Fit, 
■ SLACE139 * photoproduction 
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Figure 12. The ratio of T<i data for heavy nuclear 
targets and deuterium as measured in charged 
lepton scattering experiments in the deep inelas- 
tic region (SLAC E87 and SLAC E139) as com- 
pared to Jlab data in the resonance region. The 
scatter of the points when plotted versus x is 
large. This illustrates that the ratio does not 
scale well with the variable x. Also shown is f(x), 
which is a fit to all charged lepton scattering data 
as a function of x 



Figure 13. The ratio of Ti data for heavy nuclear 
targets and deuterium as measured in charged 
lepton scattering experiments in the deep in- 
elastic region (SLAC E87 and SLAC E139) as 
compared to Jlab data in the resonance region 
as a function of The scatter of the points 
when plotted versus £ w is smaller than when plot- 
ted versus x . This illustrates that the ratio 
scales better in the variable £ w . Also shown is 
J^updatedi^w) > which is a revised fit to charged lep- 
ton scattering data as a function of £ w . 



cross sections on free nucleons we use the follow- 
ing function obtained from a fit to SLAC data 
on the nuclear dependence of electron scattering 
cross sections j!8j . 

f{x) = 0.985 x (1 + 0.422:e - 2.745a; 2 

+ 7.570a; 3 - 10.335a; 4 + 5.422a; 5 ). (52) 

This correction shown in Fig. [14] is only valid in 
the 0.05 < x < 0.75 region. 

In neutrino scattering, we assume that the nu- 
clear correction factor for T-z, xJ-3 and 2a;J r i are 
the same. This is a source of systematic error be- 
cause the nuclear shadowing corrections at low x 
can be different for the vector and axial structure 
functions. This difference can be accounted for 



by assuming a specific theoretical model[40 a . 

14. d/u correction 

The d/u correction for the GRV98 LO PDFs is 
obtained from the NMC data for /T^ ■ Here, 
Eq. [52] is used to remove nuclear binding effects 
in the NMC deuterium T2 data. The correction 
term, S(d/u) is obtained by keeping the total va- 
lence and sea quarks the same. 

S(d/u)(x) = -0.00817 + 0.0506a; + 0.0798a; 2 , (53) 

where the corrected d/u ratio is (d/u)' = (d/u) + 
S(d/u). Thus, the modified u and d valence dis- 
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D/(n+p) (versus x) 

SLAC nuclear density Fit 




Figure 14. The total correction for nuclear ef- 
fects (binding and Fermi motion) in the deuteron, 
J 7 i/J 7 2 +P , as a function of x, extracted from fits 
to the nuclear dependence of SLAC T% electron 
scattering data. This correction is only valid in 
the 0.05 < x < 0.75 region. 



tributions are given by 

U„ 



d' = 



d v + u v S(d/u) 
1 + 



(54) 
(55) 



The same formalism is applied to the modified u 
and d sea distributions. Accidentaly, the modified 
u and d sea distributions (based on NMC data) 
agree with the NUSEA data in the range of x be- 
tween 0.1 and 0.4. Thus, we find that corrections 
to u and d sea distributions are not necessary. 

15. The Adler sum rule 

The Adler sum rules are derived from cur- 
rent algebra and are therefore valid at all val- 
ues of Q 2 . The equations below are only for the 
strangeness conserving (sc) part of the structure 
functions. These are related to the PDFs by a 
factor of cos 2 6 c . 

The Adler sum rules for the vector part of the 



structure function W% vector is given by: 



\F V (Q 2 )\ 2 + W%-™?° r (v,Q 2 )dp 
W^ tor (v,Q 2 )dv = l 



(56) 



Where the limits of the integrals are from pion 
threshold vq where W — M v + Mp to v = 00. 
At Q 2 = 0, the inelastic part of w^ vector goes 
to zero, and the sum rule is saturated by the 
quasielastic contribution \F V (Q 2 )\ 2 . Here t = 
Q 2 /4M 2 , and 

pi^^Gim+rGliiQ 2 ) 



1+T 

In the dipole approximation (discussed in an Ap- 
pendix) we have 

G V E (Q 2 ) = G D {Q 2 ) 

G V M {Q 2 ) = 4.706 G D (Q 2 ) 

G D = 1/(1 + Q 2 /M 2 ) 2 

Where M v = 0.71 (GeV/c) 2 

The Adler sum rule for y v £- axtal [ s given by: 
p 00 

\F A (Q 2 )\ 2 + W^7c a \v,Q 2 )dv 

poo 

- / W£-« w (i/,Q 9 )di/ = l 

where in the dipole approximation 

T A = -1.267/(1 + Q 2 /M 2 A ) 2 

and M A = 1.014 GeV/c 2 from reference [52], 
The Adler sum rule for w^ vector is given by: 

r\G v M (Q 2 )\ 2 + / W x V ecto >,Q 2 )^ 
W?- vector (v,Q 2 )dv=l (57) 
The Adler sum rule for yy»- axial [ s given by: 

/>00 

(1 + r)\F A {Q 2 )\ 2 + / W™ a V, Q 2 )dv 
W»-«:: al (v,Q 2 )dv = l (58) 
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The Adler sum rule for W3 is given by: 

/>oo 

2F A {Q 2 )G V M (Q 2 ) + W% n _ ac {v,Q 2 )dv 
W% p _ sc (v,Q 2 )dv = (59) 



16. Using the Adler sum rule to constrain 
the parameters of the model 



We use the Alder sum rule for w^ vector 



to 



constrain the form of the K^ff;° r ce (Q 2 ) factor for 
At low Q 2 we approximate 



•y^y— vector 



{G V E ?{Q 2 



r(G v M (Q 2 )) 2 



(1 + r) 

by G 2 D {Q 2 ) 1 and use the following K factors for 

-ys^u— vector 

KTeZeiQ 2 ) = [l-G 2 D (Q 2 )] 

where the additional parameters are obtained 
from a fit to the data. With this K™tZ.„{Q 2 ) 



factor, the Adler sum rule for WC 



L valence 
-vector 



is then 



approximately satisfied. At Q 2 = 0, the inelastic 
part of yyv- vector g 0es to zero, and the sum rule 
is saturated by the quasielastic contribution. 

In the limit of Q 2 = the Adler sum rules for 
yyaxiai anc j -yyaxmi are same. Therefore, we 

make the assumption that: 



W° B =0) = W 2 ° I! "> 



0) 



Our model is supposed to be used to esti- 
mate the cross section for W > 1.9 GeV/c 2 . 
For W < 1.9 GeV/c 2 , it is best to use a de- 
tailed resonance production model 34 35j, which 
should also includes a non-resonant background. 
However, since our model is constructed to pro- 
vide a reasonable average of the cross section 
in the resonance region (above pion threshold, 
W = 1.1 GeV/c 2 ), the transition between a res- 
onance production model and our model should 
be smooth. Note that in addition, the quasielastic 
contribution needs to be added separately. 
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Figure 15. CCFR differential cross sections 
d 2 a /dxdy on iron (in units of 10~ 38 cm 2 ) for neu- 
trinos (left) and antineutrinos (right) at an en- 
ergy of 55 GeV. The black lines in both (top and 
bottom) figures uses the default GRV98 PDFs 
without our modifications, The red line on the 
top figure is the lower limit of our predictions 
(Kaxiai = K vector ), and the blue line in the bot- 
tom figure is our prediction with the suggested 
axial K factors. 
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17. Comparison to 
Heavy Targets 



Neutrino Data on 



Figure [15] shows comparisons of the model to 
charged-current neutrino and antineutrino differ- 
ential cross section data on iron from CCFR jTSJ 
[47] at an energy of 55 GeV. On the left side we 
show the comparison for neutrino cross sections 
and on the right side we show the prediction for 
antineutrinos. In top figure the black curve is the 
prediction from GRV98 PDFs without our mod- 
ifications, and the red line is the prediction from 
our model under the assumption that the axial 
K factors are the same as the vector K factors. 
We consider this case to be a lower limit on the 
error band of our prediction. In the bottom fig- 
ure, the black curve is the also prediction from 
GRV98 PDFs without our modifications, and the 
blue line is the prediction from our model with 
our suggested axial K factors. 

Figure [16] shows the same comparisons of our 
model to charged-current neutrino and antineu- 
trino differential cross section data on iron from 
CCFR pl47] at an energy of 35 GeV. The 
CCFR data are corrected for radiative correc- 
tion. Therefore, no radiative corrections are ap- 
plied to the model predictions. However, the data 
are not corrected for non-isoscalar target (neutron 
excess). Therefore, our prediction are calculated 
for a non-isoscalar target. In both Figure [T5] and 
Figure [16] we applied the nuclear corrections from 
our fits to charged lepton scattering data to the 
model. 

We assume that the nuclear correction factor 
for J-" 2 , xJ- 3 and 1xT\. This is a source of sys- 
tematic error because the nuclear shadowing cor- 
rections at low x can be different for the vector 
and axial structure functions. This difference can 
be accounted for by assuming a specific theoreti- 
cal model [40. 

We find that our predictions are in good agree- 
ment with the CCFR |18I47| data, neutrino and 
antineutrino differential cross sections, except for 
a small difference at the lowest values of x. This 
is the region which is most sensitive to the axial 
K factors for the sea quarks, and nuclear shad- 
owing (which may in general be different neutrino 
and charged lepton scattering. 
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Figure 16. CCFR differential cross sections 
d 2 cr/dxdy on iron (in units of 10 -38 cm 2 ) for neu- 
trinos (left) and antineutrinos (right) at an en- 
ergy of 35 GeV. The black lines in both (top and 
bottom) figures uses the default GRV98 PDFs 
without our modifications, The red line on the 
top figure is the lower limit of our predictions 
{K axial — K vector), and the blue line in the bot- 
tom figure is our prediction with the suggested 
axial K factors. 
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Note however, that the lowest x region does 
not contribute much to low energy neutrino cross 
sections Monte Carlo where our model has been 
primarily used |13ll4H5j Therefore, we plan to 
use include neutrino data at low energies from 
the MINER vA [9] experiment when these data be- 
comes available in the future. 

The data are also consistent with measure- 
ments f rom CDH SW 05] and CHORUS 49j . Fig- 
ures 17 18|19|20 and 21 show comparisons of our 
model to charged-current neutrino and antineu- 
trino differential cross section data from CHO- 
RUS UK at energies of 15, 25, 35, 45, 55, 70, 90, 
110 and 130 GeV. The CHORUS data are on 
a lead target, The CHORUS data are corrected 
for radiative correction and non-isoscalar target 
(using the correction values provided in the CHO- 
RUS publication) . However, the data are not cor- 
rected for nuclear corrections for the lead target. 
In the comparison, our prediction are for for an 
iron target. Therefore, some small differences be- 
tween our predictions and the CHORUS data are 
expected because of differences in the nuclear cor- 
rections between iron and lead. 

18. Systematic errors in the application of 
the model 

The model represents all charged lepton scat- 
tering data in the deep inelastic region (W > 
1.9 GeV/c 2 ) on hydrogen and deuterium down 
to Q 2 = 0. Therefore, under the assumption 
of CVC, the vector part of the cross section in 
neutrino scattering is modeled at the 5% level. 
Therefore, the systematic error in the model pre- 
diction is dominated by other sources as described 
below. 
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• Radiative corrections: The model predicts 
neutrino cross sections at the Born level. 
Therefore, radiative corrections must be ap- 
plied to the model if it is compared to non- 
radiatively corrected neutrino or charged 
lepton data. 

• Charm sea: Since the GRBV98 PDFs 
do not include a charm sea, the charm 
sea contribution must be added separately. 
This can be implemented cither by using 



Figure 17. CHORUS differential cross sections 
d 2 a/dxdy on lead (in units of 10~ 38 cm 2 ) for 
neutrinos (left) and antineutrinos (right) at en- 
ergies of 15 and 25 GeV. The black lines in both 
(top and bottom) figures uses the default GRV98 
PDFs without our modifications, the blue line in 
the bottom figure is our prediction with the sug- 
gested axial K factors (using nuclear corrections 
for iron). 
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Figure 18. CHORUS differential cross sections 
d 2 a/dxdy on lead (in units of 10 -38 cm 2 ) for 
neutrinos (left) and antineutrinos (right) at en- 
ergies of 35 and 45 GeV. The black lines in both 
(top and bottom) figures uses the default GRV98 
PDFs without our modifications, the blue line in 
the bottom figure is our prediction with the sug- 
gested axial K factors (using nuclear corrections 
for iron). 



Figure 19. CHORUS differential cross sections 
d 2 a/dxdy on lead (in units of 10 -38 cm 2 ) for 
neutrinos (left) and antineutrinos (right) at en- 
ergies of 55 and 70 GeV. The black lines in both 
(top and bottom) figures uses the default GRV98 
PDFs without our modifications, the blue line in 
the bottom figure is our prediction with the sug- 
gested axial K factors (using nuclear corrections 
for iron). 
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Figure 20. CHORUS differential cross sections 
d 2 a/dxdy on lead (in units of 10 _3S cm 2 ) for 
neutrinos (left) and antineutrinos (right) at en- 
ergies of 90 and 110 GeV. The black lines in both 
(top and bottom) figures uses the default GRV98 
PDFs without our modifications, the blue line in 
the bottom figure is our prediction with the sug- 
gested axial K factors (using nuclear corrections 
for iron). 



Figure 21. CHORUS differential cross sections 
d 2 a/dxdy on lead (in units of 10~ 38 cm 2 ) for 
neutrinos (left) and antineutrinos (right) at en- 
ergies of 90 and 110 GeV. The black lines in both 
(top and bottom) figures uses the default GRV98 
PDFs without our modifications, the blue line in 
the bottom figure is our prediction with the sug- 
gested axial K factors (using nuclear corrections 
for iron). 



a boson-gluon fusion model, or by incor- 
porating a charm sea from another set of 
PDFs. Therefore, unless the charm sea con- 
tribution is added, the model will under- 
estimate the cross section at high neutrino 
energies and large values of v. We mod- 
eled the contribution of the charm sea us- 
ing a photon-gluon fusion model when we 
compared our predictions to photproduc- 
tion data at HERA. Note that at low neu- 
trino energies the charm sea correction is 
very small. 

Axial K factors for sea quarks: We assume 
P = 0.6 ± 0.3. With this assumption, the 
suggested axial K factors for the up, down 
and strange sea quarks and antiquarks lie 
between the assumption that the axial K 
factors for the sea are the same as the vector 
K factors for the sea (which yields a lower 
limit) and the assumption that the axial K 
factors for the sea are equal to 1 (which 
yields an upper limit) . These lower and up- 
per limits can be taken as an estimate of 
the systematic error originating from this 
source. 

Axial K factors for valence quarks: For 
the sea quarks, the vector K factor is ap- 
proximately 0.5 at Q 2 = 0.4 (GeV/c) 2 . 
For valence quarks, the vector K factor 
is approximately 0.5 at a lower value of 
Q 2 = 0.18 (GeV/c) 2 . We have made the 
assumption that the axial K factors for va- 
lence quarks is equal to 1. This assump- 
tion gives the correct neutron-proton aver- 
age cross section for the A (1232 ) reso- 
nance at Q 2 — 0. Here again, taking half 
of the difference between assuming that the 
valence axial K factor is equal to 1, and the 
assumption that the valence axial K factor 
is the same as the valence vector K factor 
can be taken as the systematic error from 
this source. 

Nuclear corrections: We have assumed that 
the nuclear corrections for the axial and 
vector case are the same, and equal to 
the nuclear corrections measured in charged 



lepton scattering. We have also assumed 
that they are the same for all three struc- 
ture functions, and that these corrections 
can be applied by using a simple fit in 
x w . The systematic error from this source 
can be reduced by using specific theoreti- 
cal models [40] to account for the differences 
in the nuclear corrections between neutrino 
and charged lepton scattering as a function 
of Q 2 and x for various nuclear targets. Be- 
cause of the uncertainty in the nuclear cor- 
rections, we chose to model neutrino and 
charged lepton scattering on hydrogen and 
deuterium only, and we did not include any 
neutrino data (which are mostly on nuclear 
targets) in our fits at this stage. 

• A better determination of the axial K fac- 
tors for the sea quarks requires including 
neutrino data (e.g. CCFR, CDHS, CHO- 
RUS, NOMAD and MINERvA) on nuclear 
targets in our fit. The nuclear corrections 
at small x for vector and axial structure 
functions are different in general. There- 
fore, it is not possible to separately deter- 
mining the axial K factor without assum- 
ing a specific model for the neutrino nu- 
clear corrections [40 . This study will be pre- 
sented in a future communication. 

19. Using the model near Q 2 = 

Near Q 2 = and x = there are terms in 
the expression for the cross section that are pro- 
portional to the square of the muon mass as 
In the low Q 2 region it 



shown in equation 35 



is more convenient to use the structure functions 
T\ = MWi, T = vW 2 , and J" 3 = vW 3 , instead 
of 2xJ r i and xT 3 which are more useful at higher 



The expression for the scaling variable 



2x{Q 2 + Mt 2 + B) 



Q 2 [l + + AM 2 x 2 /Q 2 ] + 2Ax ' 
can be written alternatively as 
(Q 2 + M f 2 + B) 



Mv[l + ^Jl + Q 2 /v 2 ]+A' 



(61) 



(62) 
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The structure function 2xT\ should be divided 
by 2x to yield T\ and the structure function xJ-'a 
should be divided by x to get J 7 ^. Note the T\ 
and J-~3 (or alternatively Wi and W3 ) structure 
functions do not diverge at Q 2 = 0. This because 
2xJ r i and xJ-3 are proportional to Q 2 near Q 2 = 
(because of the various K factors) . 

For example, the valence d quark component of 
J-3 is given by: 



J". 



■valence— down 



(Q\w) 



2 /Q 2 [l-Gi(Q 2 )] 



At low Q 2 the vector K factor for valence 
quarks [1 — G 2 D (Q 2 )] can be written as Q 2 / (Q 2 + 
0.178), which yields. 



T. 



■valence— down / s~\2 



(Q 2 ,W) = 



2Mi- 



V(Q 2 +CxF3)(Q 2 +0.178) 

xK Lw t 11 , (f o 2 ) ( g 2 ±sn!!^ 

20. Using the model to determine the flux 
of low energy neutrino beams 

It is very difficult to make a good measurement 
of the neutrino and antineutrino flux as a function 
of energy in a broad band beam such as the NUMI 
beam. The following technique can be used to 
extract the flux directly from the data. 

The model gives a very good prediction of 
the neutrino differential cross sections in a re- 
stricted region of phase space. This region cor- 
responds to W > 1.9 GeV/c 2 (the inelastic re- 
gion) Q 2 > 1 (GeV/c) 2 (where vector and axial 
structure functions are close to each other) and 
x > 0.1 (where nuclear effects in neutrino and 
charged lepton scattering should be very similar) . 
The data in this restricted region can be used to 
determine the neutrino and antineutrino flux by 
comparing it to our model predictions. 

At very low energies there may not be many 
events in the region of validity of our model. How- 
ever, we can cross calibrate between high and low 
energies by using events at a fixed range of W and 
Q 2 < 1 (GeV/c) 2 . This is because as can inferred 
from equation [35] the differential cross section for 



a fixed range of W and low Q 2 is dominated by 
W2 and is relatively independent of neutrino en- 
ergy, and the small energy dependent corrections 
given in equation [35] can be calculated. We will 
present a study of these two techniques in a future 
communication. 

21. Appendix - Quasi-elastic Neutrino- 
Nucleon Scattering 

In order to calculate the total neutrino and an- 
tineutrino total cross sections, we need to add the 
quasi-elastic cross sections. Here we use the no- 
tation of C.H. Llewellyn Smith [51] 

q 2 = ?o -4 = -4£o£' sin 2 l = -Q 2 ■ 

The hadronic current for QE neutrino scatter- 
ing is given by [51] 

<p(Pi)\Jt l n (Pi) >= 



jxH(q 2 )~ 

7a75-7 7 a(<? 2 ) 



2NI 



u(pi 



where q — k v — k^, £ = (fi p — 1) — (i n , and M = 
(m p + m n )/2. Here, \i v and \x n are the proton 
and neutron magnetic moments. We assume that 
there are no second class currents, so the scalar 
form factor JFy and the tensor form factor T\ 
need not be included. Using the above current, 
the cross section is 



da v 



M 2 G 2 F cos 2 e c 



dq 2 



A(q 2 ) T 
where 
A(q 2 



8ttE 2 

( S -u)B(q 2 ) , C{q 2 ){ S -uf 



M 



2 2 
rrr — q z 



AM 2 
M 2 



M 2 



M 



\n\ 2 



M 



4M 



Aq 2 ReT v *jT 2 
M 2 



B{q 2 ) 



M 2 



ReF* A {F v +(iF v ), 



2G 



C(q 2 



\rt\ 2 



q 

M 2 



in 



Although we have have not shown terms of or- 
der (mi/M) 2 , and terms including Fp(q 2 ) (which 
is multiplied by (m;/M) 2 ), these terms arc in- 
cluded in our calculations [51].) The form factors 
Fy{q 2 ) and £,F 2 (q 2 ) are given by: 



iF 2 { q 2 



G V E (q 2 ) 



AM 2 



G\M 2 )-G V E {q 2 ) 



AM 2 



From conserved vector current (CVC) G^(q 2 ) 
and G\[{q 2 ) are related to the electron scatter- 
ing form factors G p E (q 2 ), G E (q 2 ), G p M (q 2 ), and 

G n M (q 2 y- 

G E (q 2 ) = G E (q 2 )-G E (q 2 ), 

G v M { q 2 ) = G p M {q 2 )-G n M {q 2 )- 

At low Q 2 the axial form factor Ta can be ap- 
proximated by the dipole form 



?A{q 2 ) 



1 



.9,1 



Mi 



In our analysis we apply BBBA2008 
corrections [52] to the the dipole parametriza- 
tion of the axial form factor as described in 
reference [S5]. 

The pseudoscalar form factor Tp is related to 
Fa by PC AC and is given by: 



2M 2 T A (q 2 ) 
Ml -a 2 ■ 



In the expression for the cross section, Tp{q 2 ) is 
multiplied by (mi/M) 2 . Therefore, in muon neu- 
trino interactions, this effect is very small except 
at very low energy, below 0.2 GeV. 

At low Q 2 a reasonable description of the pro- 
ton and neutron elastic form factors is given by 
the dipole approximation Gd — l/(l + <5 2 /0.71) 2 . 



The dipole approximation is a lowest-order at- 
tempt to incorporate the non-zero size of the pro- 
ton into the form factors. If we also assume that 
the proton has a simple exponential spatial charge 
distribution: 



p(r) 



= Poe- r/ro 



(63) 



where ro is the scale of the proton radius. The 
form factors are related, in the non-relativistic 
limit, to the Fourier transform of the charge and 
magnetic moment distribution. If it is also as- 
sumed that the magnetic moment distribution 
has the same spatial dependence as the charge 
distribution (i.e., form factor scaling), we get the 
dipole approximation to the form factors: 

G D (Q 2 ) ' 



(i + Q 



2„2\2 



(l + QVO-71) 2 

G p E (Q 2 ) 

G p M {Q 2 )/» P 

G n M (Q 2 )/^n 



and G n E = 



0. Here fi p = 2.7928, and fj, n = 
1.913. Previous measurements of e-p and e-n elas- 
tic scattering have indicated a best fit value of 
r 2 = (0.24 fm) 2 = 1/0.71 (GeV/c) 2 , indicating an 
rms radius of y/ (r 2 ) « 0.81 fm. Measurements of 
G P E and G P M agree with the dipole approximation 
only for Q 2 < 1 (GeV/c) 2 . The dipole approxi- 
mation yields. 

G V M {Q 2 ) = 4.706 GD(Q 2 ) 

Using r = Q 2 /M 2 we also get 



fUQ 2 



G E (Q 2 



tG\j 



(Q 2 



l 



In our calculation of the quasi-elastic cross 
sections we apply BBBA2008 corrections to the 
dipole parametrization of the vector form factors 

(S2j. 

Iii evaluating the quasielastic contributions to 
the Adler sum rules we use: 



-\ j tv— vector 
ln—elasti 



s<r - ^r)G v M {Q 2 ) 2 
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W 



in-l a L = S(y + r)\T A (Q 2 ){< 



if 



1/ — vectov 
In— elastic 

-\jtTV— axial 
27i— elastic 



2 Ml5(u-$-)\FUQ 
v v 2M J] VW 



2\|2 



Mt 

V 



22. Appendix -Results with GRV94 PDFs 

and x w 

For completeness we describe our earlier analy- 
sis (1011 1] in which we used another modified scal- 
ing variable [H] x w with GRV94 PDFs(instead 
of GRV98) and simplified K factors. In that anal- 
ysis we modified the leading order GRV94 PDFs 
as follows: 

1. We increased the d/u ratio at high x as de- 
scribed in appendix A [16) . 

2. Instead of the scaling variable x we used the 
scaling variable x w = (Q 2 + B)/(2Mv + A) 
(or =x{Q 2 + B)/{Q 2 + Ax)). This mod- 
ification was used in early fits to SLAC 
data [57]. The parameter A provides for 
an approximate way to include both target 
mass and higher twist effects at high x, and 
the parameter B allows the fit to be used all 
the way down to the photoproduction limit 

(0 2 =o). 

3. In addition as was done in earlier non-QCD 
based fits |24I25] to low energy data, we 
multiplied all PDFs by a factor K=Q 2 / 
(Q 2 +C). This was done in order for the 
fits to describe low Q 2 data in the photo- 
production limit, where T2 is related to the 
photoproduction cross section. 

4. Finally, we froze the evolution of the GRV94 
PDFs at a value of Q 2 = 0.24 (for Q 2 < 
0.24), because GRV94 PDFs are only valid 
down to Q 2 = 0.23 (GeV/c) 2 . 

As was done for GRV98, in the GRV94 analysis, 
the measured structure functions were also cor- 
rected for the BCDMS systematic error shift 20J 



and for the relative normalizations between the 
SLAC, BCDMS and NMC data [THIT7| . The deu- 
terium data were corrected for nuclear binding 
effects |16ll7j . A simultaneous fit to both pro- 
ton and deuteron SLAC, NMC and BCDMS data 
(with x > 0.07 only) yields A=1.735, B=0.624 
and C=0.188 (GeV/c) 2 ) with GRV94 LO PDFs 
(X 2 = 1351/958 DOF). Note that for x w the 
parameter A accounts for both target mass and 
higher twist effects. 

In our studies with GRV94 PDFs we used the 
earlier R wor id fit [19] for R ncp and R cp . R wor id is 
parameterized by: 



Rworld{x, Q 2 > 0.35) 



0.0635 



x,Q 2 



log(Q 2 /QM) 
0.5747 0.3534 



Q 2 



0.09' 

(64) 



where 9 = 1. 



12Q 2 



0.125^ 



Q 2 + 1.0 ~ 0.125 2 +ic 2 / Tfl e R W orld 

function provided a good description of the 
world's data for R at that time in the Q 2 > 0.35 
(GeV/c) 2 and x > 0.05 region (where most of 
the R data are available). However, for elec- 
tron and muon scattering and for the vector part 
of neutrino scattering the R wor id function breaks 
down below Q 2 = 0.35 (GeV/c) 2 . Therefore, we 
freeze the function at Q 2 = 0.35 (Gev/c) 2 . For 
electron and muon scattering and for the vector 
part of T\ we introduce a K factor for R in the 
Q 2 < 0.35 (Gev/c) 2 region. The new function 
provides a smooth transition for the vector R 
(we use R V ector—R e /fi) from Q 2 = 0.35 (Gev/c) 2 
down to Q 2 = by forcing R vec tor to approach 
zero at Q 2 — as expected in the photopro- 
duction limit (while keeping a 1/Q 2 behavior at 
large Q 2 and matching to R wor id at Q 2 = 0.35 
(Gev/c) 2 ). 



R ve ctor(x,Q 2 < 0.35) = 3.207 x 



Q 2 



Q 2 = 0.35). 
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